Priedel oscillations in a two-band Hubbard model for CuO chains 
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Friedel oscillations induced by open boundary conditions in a two-band Hubbard model for 
CuO chains are numerically studied. We find that for physically realistic parameters and close 
to quarter filling, these oscillations have a 2/cf modulation according with experimental results on 
YBa2Cua07_{. In addition, we predict that, for the same parameters, as hole doping is reduced 
from quarter filling to half filling, Friedel oscillations would acquire a 4/cf modulation, typical of 
a strongly correlated electrons regime. The 4&f modulation dominates also in the electron doped 
region. The range of parameters varied is very broad, and hence the results reported could apply 
to other cuprates and other strongly correlated compounds with quasi-one dimensional structures. 
On a more theoretical side, we stress the fact that the copper and oxygen subsystems should be 
described by two different Luttinger liquid exponents. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.10.Pm, 71.10.-w 



I. INTRODUCTION 

The application of scanning tunneling spectroscopy 
(STM) techniques to the CuO chain planes in 
YBa2Cu307-«5 (YBCO) has provided new insights on 
the physics of these quasi-one dimensional (ID) elec- 
tronic structures. These experimental studies have shown 
in real space the presence of charge modulations along 
the chains at low temperature, inside the supercon- 
ducting phasei^ These charge oscillations had been in- 
ferred earlier from neutron scattering experiments on 
the same compound. 3 Taking together results from both 
band calculations^ and angle-resolved photoemission 5 
(ARPES) indicating a chain filling close to one quarter, 
it turns out that the observed charge modulations have 
a 2k p wavenumber. 

Although earlier STM studies^ had given indications 
of the presence of charge oscillations on the chains, the 
present interest on these features comes from a more re- 
cent set of experiments exploring the interplay between 
the chain plane and the Cu02 plane in YBCO. On one 
side, the chains (running along the b direction) may 
contribute to the in-plane resistivity anisotropy in this 
compound^ It is also important to notice that stripes in 
Cu02 planes also run along the b direction, may be not 
coincidentally. On the other side, there have been exper- 
imental indications of a superfluid density induced on the 
chains due to the proximity to the Cu02 planes&S Al- 
though the interpretation of this last set of experiments 
is still controversial 5 ' 10 , one recent theoretical studj>i» 
offers an explanation of STM results! based on this sce- 
nario of proximity induced chain superconductivity (SC). 
In this model, the resonances in the chains arise from the 
interference between magnetic impurities in the chains. 

Alternatively, we believe it is necessary to search for 
explanations of these experimental results based on mod- 



els which capture the intrinsic electronic interactions on 
the chains. Eventually, in this kind of models, the SC on 
the Cu02 planes could modify at an effective level the 
coupling constants of the intrachain interactions. This 
approach connects the present study to the more general 
theoretical interest in various types of charge inhomo- 
geneities appearing in strongly correlated electron sys- 
tems in low spatial dimensions. In fact, in addition to 
the relevance of the presence of charge oscillations in the 
chains to the physics of the Cu02 planes, it has been 
emphasized that precisely the assumed strong Coulomb 
interactions in the cuprates would naively favour a AUf 
instead of the observed 2k p wavenumber of the charge 
modulations^ 

Since CuO chains are cut by oxygen depletion, we 
will try to describe the charge modulations as general- 
ized Friedel oscillations (FO) starting at the open end 
of the resulting fragments. A very recent numerical and 
analytical studyi^ on the ID one-band Hubbard model 
has indeed found a crossover from 2kp to Mzf FO as the 
strength of the onsite Coulomb interaction U is increased. 
We are going to examine a two-band Hubbard model, ap- 
propriate for the CuO chains, where the situation is more 
complex due to different Coulomb repulsions on Cu and 
O ions, in addition to the charge transfer between them^ 
It is interesting to notice that ARPES shows^ a chain's 
Fermi surface in agreement with band calculations and 
at the same time the dispersion along the chain direction 
agrees with the holon band predicted by strongly corre- 
lated theories.— Hence, although the Coulomb repulsion 
on Cu ions is large, it does not automatically imply a 
strongly correlated behavior on any given physical quan- 
tity. In particular, we will show a complex dependence 
of 2k f and 4kp FO with respect to the Coulomb inter- 
actions and charge transfer parameter. 

The model here studied does not include electron- 
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lattice coupling which is another possibility of inducing 
2kp instabilities. Experimentally^, it has been suggested 
that standard lattice phonons are unlikely to reproduce 
the observed features. Finally, both in order to help the 
understanding of our results and in order to extend the 
scope of this study, we will consider a somewhat large 
variation of interaction parameters, and in addition to 
quarter-filling, dopings close to half-filling, and also elec- 
tron doped chains will also be studied. Besides, we will 
consider open ends occupied with O ions and with Cu 
ions. 

This paper is structured as follows. In Section [H] we 
describe the model and method of calculation, and we 
show the numerical results for density oscillations. The 
theoretical and experimental implications of these results 
are discussed in Section UTTl 



II. 2k F /Ak F PHASE BOUNDARY 
A. Model and method of calculation 

The model studied is the ID two-band Hubbard model 
defined as: 

H= - tpd^icl^^tv + h.c.) 

+ Ud ^ »2i-iT ra 2i-ij, + Up ^ ri 2 i-\n 2 ii 

i i 

n 2 i-in 2 i + A^n 2i (1) 

i i 

where, cj CT creates a hole with spin a at site j, nj a — 

C ]a C j<? 



, n.j — rij^+rij^. The sums run over i = 1, . . . , L, L: 
number of unit cells, and Cu (O) ions occupy odd (even) 
sites. t p d = l- All our results shown below correspond to 
Ud = 8, Upd = 0, and U p , A < Ud- The number of holes 
(Nh) is larger than L for hole doped systems in which case 
the filling is computed as nu = (2L — Nh)/L, or smaller 
than L for electron doped systems in which case the filling 
is computed as n e = Nh/L. At half filling = n e = I. 
The numerical technique employed in the present study 
is the density matrix renormalization group algorithm 
(DMRG)ii^ The most important quantity measured by 
this technique is the on-site charge density, n(i) — (rii). 
Most results were obtained with a truncation number 
in = 200. For several cases we examined also m = 300 
and m = 400 with almost indistinguishable results. 

In principle, we are interested in the effect of oxygen 
depletion which implies that CuO chains are cut at oxy- 
gen sites. These depleted oxygens could be modeled by 
just imposing open boundary conditions on chains with 
odd number of sites with the two end sites correspond- 
ing to Cu-ions. However, it is technically convenient to 
use chains with even number of sites which implies that 
one end is occupied by an oxygen ion and the other end 
site by a copper ion. The open end with an O ion would 



correspond to a nonmagnetic impurity, for example a Zn 
ion replacing a Cu ion>i& If one assumes that the Friedel 
oscillations from both ends do not have nonlinear super- 
position, i.e. both oscillation are independent of each 
other, then the use of even number of sites allows us the 
simultaneous study of both types of impurities in a single 
chain. 

Since, as said in the Introduction, Friedel oscillations 
appear due to oxygen depleted sites acting as impurity 
centers in an otherwise metallic chain, then they will fol- 
low a power law behavior predicted by Luttinger liquid 
theory?i£ 



n(r) — tiq ^ a\ 



cos(2kpr + cf)i) 

r (l + K p )/2 



cos(4fcF?- + <fa) ,_, 
° 2 J2K ' ( 2 ) 



with r the distance from the impurity. This expression 
can be derived from the asymptotic form of the charge 
correlations and K p is the interaction dependent corre- 
lation exponent^ Numerical studies have verified that 
FO follow this form for the one-band Hubbard and t-J 
models on chains and ladders 12 , and for the Kondo lattice 
model on chains^ In the following, we will fit our data 
using Eq. . The complete fitting procedure is detailed 
in the Appendix. The key point of this procedure is that 
we simultaneously fit the FO starting from the left chain 
end (Cu site) and the one which starts from the right 
chain end (O site). 



B. Quarter filling 

We start by examining the electron density n{i) as a 
function of the position obtained by DMRG for L = 64, 
Nh = 96, which corresponds to quarter filling, = 0.5, 
as a function of U p and A. At this filling, the density 
modulation has period of four unit cell spacings for the 
2kF component. This period can be clearly seen in the 
Friedel oscillations shown in Fig. ^ (U p = 2) in both 
Cu and O ions, at A = 2. These density oscillations 
extend appreciably over the whole chain. For A = 5, the 
density oscillations on both type of ions present a period 
of two which corresponds to a 4fcp wave. Also shown are 
the oscillations for A = 4, close to the crossover between 
both regimes. 

By applying the fitting procedure previously discussed, 
we are able to obtain the mean amplitudes for the 2kp 
and 4k f components of the Friedel oscillations on Cu 
or on O ions starting at site 1 (occupied by a Cu ion), 
corresponding to an impurity on an O ion, and the one 
starting at site 64 (occupied by an O ion), correspond- 
ing to an impurity on a Cu ion. These mean amplitudes 
are defined by Eqs. I|A4(1 and (|A5I) . The inclusion of 
the density oscillation close to the ends introduces some- 
what large numerical errors so we discarded the first three 
points near the ends. The results are depicted in Fig. |2] 
The error bars are about the size of the symbols used. 
By discarding different number of initial sites of course 
the amplitudes vary. However, their crossing point does 
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FIG. 1: Density at site i for L = 64, n h = 0.5, U d = 8,U P = 2, 
U p d = and A = 2, 4 and 5 from bottom to top. O sites: 
open symbols, Cu sites: filled symbols. The curves have been 
arbitrarily shifted for the sake of clarity. 
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FIG. 2: A dependence of the mean amplitude of 2&f (open 
symbols) and 4fc_F (filled symbols) components for Ud = 8, 
L = 64 and various values of J/ p at ri|, = 0.5. (a) O density os- 
cillations, O-impurity, (b) O density oscillations, Cu-impurity, 
(c) Cu density oscillations, O-impurity, (d) Cu density oscil- 
lations, Cu-impurity, 



FIG. 3: Phase boundary between the 2Uf and 4fcp Friedel 
oscillations at quarter filling for Ud = 8, L = 64, (a) O- 
impurity; (b) Cu impurity. Open circles: boundary of Cu 
oscillations, filled circles: boundary for O oscillations. Dotted 
lines are lines of constant J/t e ff of the effective one band t-J 
model (see text). The gray circles correspond approximately 
to the physical region of CuO chains. 



rized in the "phase" diagram shown in Fig. [3J The re- 
gion of A below the crossover corresponds to the 2fcf- 
dominated regime while the region of A above it corre- 
sponds to the 4fcF-dominated regime. The first remark- 
able feature is that the crossover in the O-ions oscillations 
takes place at lower values of A than the crossover in the 
Cu-ions oscillations. 

To help the interpretation of these results, we have 
overimposed on this figure lines of constant J/t e f / of an 
effective one-band t-J model with)3I 



A 2 
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Ud 



2A + U t 



and 



f 2 
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(3) 



(4) 



not change appreciably. Various features can be rapidly 
appreciated: a larger amplitude on O sites than on Cu 
sites for the dominant component; the suppression of the 
2kp amplitude with increasing U p for O-site oscillations; 
larger (smaller) O (Cu) amplitudes for Cu-impurity than 
for O-impurity. 

The crossing of these mean amplitudes as A is in- 
creased for a given value of U p is adopted as the cross- 
ing of the 2fci?-dominated regime to the 4fc^-dominated 
regime. Other criteria could be eventually adopted but 
we found that this criterion faithfully reproduces the 
change of oscillation period which is the quantity mea- 
sured in STM experiments. The results can be summa- 



We have fixed the proportionality constant in such a way 
that the physical point for the CuO chains corresponds 
to J/t e ff = 0.35 — 0.4, taking into account some disper- 
sion on the values reported for the Coulomb interactions 
which are approximately, Ud ~ 8 (the value adopted 
in the present study), U p ~ 3 — 4, A ~ 3 (in units 
of t p d)- At a fixed hole density, each line of constant 
J/teff (indicated in Fig. 01 corresponds in the effective 
ID t-J model to a line of constant K p , which have been 
previously estimated^ Then, it can be appreciated in 
Fig. |21 that the IkpjAkp boundaries are roughly parallel 
to these lines with the largest deviation corresponding to 
the O oscillations, Cu impurity. In addition, as it can be 
seen in Fig. |3fa). for the physical region corresponding 
to the CuO chains (indicated with a gray circle), the 




FIG. 4: Density at site i for L = 64, n h = 0.875, U d = 8, FIG. 5: Same as Fig. but for n h = 0.875. 

Up — 0, U p d = and A = 1, 3 and 5 from bottom to top. 
O sites: open symbols, Cu sites: filled symbols. The curves 
have been arbitrarily shifted for the sake of clarity. 



FO generated by oxygen depletion have 2kp modula- 
tion in agreement to experimental results near quarter 
filling <* In the case of Cu impurity (Fig. Efb)), although 
the 2fcF/4fcp boundary is shifted to lower values of A, 
the FO are predicted to have also a,2kp character in the 
physical region. 



C. n h = 0.875 filling 

While in the undoped compound, YE^CusCv, the 
chain filling is = 0.5, upon oxygen depletion it in- 
creases toward half-filling, which is reached at S w 0.36 
fRef. l2lf> . Hence, in this subsection we consider the case 
of i%h — 0.875. Although increased O depletion would 
produce on average shorter CuO segments, in order to 
reduce the number of varying parameters and to facil- 
itate the comparison, we study the same chain length 
as in the previous subsection, i.e. L = 64. The filling 
rih = 0.875 then corresponds to Nh — 72 holes. In this 
case, 2Uf — 7n/8 and 4k f = 7f/4. 

Let us start by examining in real space density oscil- 
lations corresponding to this filling. They arc illustrated 
in Fig. 21 The oscillations for small A, corresponding to 
2kp can be regarded as period 2 oscillations with kinks, 
while the oscillations at large A, corresponding to Ak f 
have a well defined period 8 unit cell spacings. 

As we did in the previous subsection, we have fitted 
the density oscillations for both Cu and O ions and the 
results for the mean amplitudes of the 2k p and 4k p com- 
ponents are shown in Fig. EI By comparing with the pre- 
vious results at quarter filling, it can be concluded that 
the amplitudes of the dominant component are in gen- 
eral much larger in this case of = 0.875. Again O 
oscillations have in general larger amplitudes than the 
Cu ones. It can be also observed that the effect of U v is 



relatively small in the present case. Also there are virtu- 
ally no differences between the cases of O-impurity and 
Cu-impurity. 

Again, the crossing points between 2fcp and 4fcp com- 
ponents of the mean amplitudes lead to the "phase" dia- 
gram shown in Fig. El As anticipated by the previous fig- 
ure, the boundary between the 2kp (low A) and the 4kp 
(high A) regions has been considerably shifted to lower 
values of A with respect to the quarter filled case. As for 
nh = 0.5, these boundaries also follow approximately the 
lines of constant K p of the effective ID t-J model. How- 
ever, in contrast with that case, there is much less dif- 
ference between the boundaries for O depletion-induced 
FO (Fig. EJa)) and the boundaries for the FO induced by 
nonmagnetic substitution of Cu ions (Fig.EJb)). In any 
case, as in the quarter filled case, the most remarkable 
feature is again that the 2k f to 4kp crossover for oxygen 
FO takes place at smaller values of A than the one for 
Cu FO. Notice that for this filling the physical region of 
parameters (gray circles) falls now on the 4k f side of the 
diagram. 

Finally, by further oxygen depletion, 5 > 0.36, the 
chain would formally go into the electron doped region. 
The 2fcp/4fci? boundary between Friedel oscillations on 
Cu ions in the electron doped, n e — 0.875, system is 
also included in Fig. El This result confirms the trend 
suggested by the cases nh = 0.5 and nh = 0.875, that 
by reducing the number of holes this boundary shifts to 
lower values of A for a fixed U p . On the other hand, the 
occupancy on O sites is very low and we observed 2k f 
oscillations for all the values of parameters examined. It 
is reasonable to speculate that by further reducing the 
number of holes the presence of 4/cf oscillations on Cu 
ions will eventually disappear as well, i.e. the effects 
of strong electron correlations are expected to be more 
important close to half-filling. 
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FIG. 6: Phase boundary between the 2&f and 4fci? Friedel os- 
cillations for Ud = 8, L — 64, rih = 0.875. (a) O-impurity; (b) 
Cu impurity. Open circles: boundary of Cu oscillations, filled 
circles: boundary for O oscillations. Open triangles corre- 
spond to Cu oscillations in electron doped, n e = 0.875, chains. 
The gray circles correspond approximately to the physical re- 
gion of CuO chains. 



D. Global properties 

Let us examine some global properties like the total 
ground state energy Eq and the average amplitude of FO. 
Since this energy is monotonically increasing with A, for 
all the other parameters fixed, it is more meaningful to 
analyze the quantity obtained by subtracting from Eq 
the quantity Ano where no is the average filling on O 
sites. The average of the amplitude of the FO on Cu sites 
is defined as, 

Acu = j^^2\n{i) - n Cu \ (5) 

i 

where the sum runs over Cu sites and nc u is the average 
filling on Cu sites. To be consistent with the fitting pro- 
cedure and hence to allow comparison with Figs. [5] and 
we discarded the first three sites close to both chain 
ends. The resulting quantity is shown in Fig. EJa),(b) 
for rih = 0.5 and rih = 0.875 respectively. It can be no- 
ticed that it presents a minimum at a value of A close to 
the 2fci?/4fcf boundary of Cu FO at a given value of U p . 
As it can be seen in Fig. [7fc) , (d) , the corrected energy, 
Eq — Ano also presents a minimum at approximately the 
same values of A for — 0.5. At rih — 0.875, the min- 
imum of this quantity is located at a lower value, closer 
to the 2kp/Akp boundary of oxygen FO. 



III. DISCUSSIONS AND CONCLUSIONS 

The two most important results of the present study, 
which can be deduced from Figs. Eland HJ|are: (i) the 



FIG. 7: (a) Average amplitude of Friedel oscillations on Cu 
sites for rih = 0.5. (b) Same for nh = 0.875. (c) Corrected 
ground state energy for rih = 0.5. (d) Same for rih = 0.875. 



2k p to 4/cp crossover takes place at smaller values of A 
for given U p as doping is reduced from quarter filling to 
half filling, and (ii) for hole doped chains this crossover 
takes place at smaller A for the density oscillations on O 
ions than for the ones on Cu ions. Let us discuss these 
two features in detail. 

In the first place, the fact that the Ikpl^kF bound- 
aries quite likely correspond to lines of constant K p of 
an effective ID t-J model suggests that both the approx- 
imate mapping to this model and the whole fitting pro- 
cedure are physically correct. One should notice that 
this effective one band model is obtained essentially by 
projecting out the oxygen sites, i.e. is a model valid 
mostly for Cu ions. Besides, it becomes less valid as one 
moves away from half filling. It is not obvious then that 
also the boundaries for oxygen oscillations follow approx- 
imately the lines of constant J/t e f / although with some 
important deviations, specially at quarter filling. The 
dependence of K p with J/t and hole density for the ID 
t-J model 22 , is consistent with the shifting of the Cu 
2&f/4£;f boundary to lower values of A as hole density 
rt/j — > 1, and in turn this is consistent with the intuitive 
notion that strongly correlated electrons regimes become 
more important close to half-filling. 

However, this mapping does not explain by itself why 
in the case of hole doping, the 2kp / 4kF boundaries for O 
oscillation occur at smaller values of A for any given U p , 
i.e. why the O-subsystem enters in the strong correlation 
regime at smaller interactions than the Cu-subsystem. 
To understand this feature, it would be necessary to ana- 
lyze effective models obtained by projecting out Cu sites 
and retaining O sites. These effective models, the so- 
called "spin-fermion" or Kondo-Heisenberg models*^ are 
considerably more complicated and related available re- 
sults are limited^ For hole doped cuprates, i.e. with 
holes in excess of the number of Cu ions, it can be as- 
sumed that Cu ions are mostly single occupied, and in 
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this case our results suggest that the 2kp to <lkp crossover 
is "driven" by the oxygen oscillations. That is, the O sub- 
system is more susceptible to strong correlations than the 
Cu subsystem. On the other hand, for the case of elec- 
tron doped chains, i.e. with fewer holes than Cu ions, the 
2kp to Akp crossover takes place at small A for Cu den- 
sity oscillations while for O FO have 2k p in the range 
of parameters studied, and in this sense, the crossover 
seems "driven" by the Cu ions. Hence, the nearly empty 
O band in electron doped chains is symmetric to the half- 
filled Cu band in hole doped cuprates. 

There are two factors affecting the magnitude of the 
mean amplitudes of the 2kp and 4fcp components: their 
amplitude at the origin, and their power law exponents. 
The amplitudes at the origin (A 2 , A 3 , B 2 , B 3 , in Eqs. 
l|A"Tj) and EU) may depend on the formation of bound 
states between holes and impurities. With respect to the 
power law exponents, one could assume that the FO on 
Cu and on O sites are determined by two different K p 
since they obey different effective models in addition of 
having different hole densities. Unfortunately, although 
the fitting to the oscillations using Eqs. (|A1|> . ljA2() . (|A3|) 
are very good, the determination of K p from the fitted 
values is very noisy to allow reliable conclusions. In this 
sense, an independent study of periodic chains, with a 
careful study of finite size effects, would be necessary to 
confirm this possibility. 

Order of the crossover. It seems to be a first order in 
the sense that there are no other oscillation wavenumber 
between 2kp to The evidence for this is rather in- 

direct since our method of fitting does not make room 
for variable wavenumber, and a Fourier (or windowed 
Fourier) transform would mix the Friedel oscillations 
starting form the left and right ends of the open chain. 
The indication of a first order type of crossover comes 
from the fact that the overall average amplitude of the 
density oscillations are minimum at the crossing point 
(Figs. [J Hand EJa),(b)). This suggests that the sys- 
tem becomes "frustrated" at the crossing point due to the 
competition of 2k p and -ikp modulations being unable to 
develop a modulation at another wavenumber. The be- 
havior of the corrected ground state energy, E — Ano, 
(Fig- Etc), (d)) is consistent with this interpretation: its 
minimum is located close to the point where Ac u is also 
minimum presumably because it gains energy due to de- 
localization. Hence, this quantity also suggests that there 
are no FO with a modulation intermediate between 2k p 
and Akp. 

Predictions. The main prediction of the present study 
is that CuO chains in YBa2Cu307_,5 would undergo a 
2k p to Akp crossover in the Friedel oscillations induced 
by oxygen depletion or by nonmagnetic substitution of 
Cu ions as doping moves from quarter to half filling and 
eventually into the electron doped region. Additionally 
we predict that experiments like STM, probing O sites, 
would detect a Akp modulation while experiments such 
as neutron scattering would still see a 2k p modulation 
on Cu sites. More generally, the modulation of FO can 



be considered as a sensible tool to detect the sometimes 
subtle presence of strong electron correlations in quasi- 
1D systems. Finally, if the observed modulations in CuO 
chains are essentially Friedel oscillations and the presence 
of superconductivity on the planes plays a minor role, 
then similar modulations should be observed in the non- 
superconducting compound PrBa2Cu307 fRef. 124 
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APPENDIX A: FITTING PROCEDURE 

In general, the oscillation starting from the Cu-edge 
is different to that from the O-edge. The total oscil- 
lation, as obtained numerically, is the superposition of 
both oscillations. Assuming these oscillation follow the 
Luttinger power law expressions, 

1 cos(2fc F r + (fti) cos(4fc F r + fc) 
L(r) = A +A l -+A 2 — +A 3 — , 



and 

R{r) = B Q + B l 
+ B 3 



1 



B- 



(Al) 

cos(2kp(L — r) + cpi) 



L-r ' ~" (L-r)n 
cos(4:kp(L — r) + ip 2 ) 
(L - ryn* ' 



(A2) 



where L is the system size, then our fitting function will 
be, 



S(r) = L{r)+R{r) 



(A3) 



with no < r < L — hq. tiq is the number of sites from 
each edge that are cut out to avoid wild oscillations. Be- 
low, the constant terms, Aq and B , are neglected, since 
a uniform component is subtracted before the fitting pro- 
cess. 

The mean amplitudes of the 2kp and Akp components 
of the Friedel oscillations starting from the left and from 
the right edges of the chain are defined as, 



.4 



2k F .L 



A 



4k F ,L 





A 2 


/; 


-no ^ 

dr 

r 7l 


L 


- 2n 






A 3 


r 


-no ^ 

dr 

J- 72 


L 


- 27i 


J n 



(A4) 



7 



and 



r 


-no I 

dr 


'no 




r 


-no ^ 

dr 


J no 





= i^k i dr ~ (A5) 

respectively. 

In the fitting process, Eq. (|A3(I is rewritten as, 

Qf \ , 1 , a, cos{2k F rj) sm(2k F n) 
&(n) = Ai— + A 2 v A 2 

'« ' i 1 i 

A , cos(4k F n) „ sm(4h F ri) 1 
+ A3 To rA* ^ V Bi- 



cos(2fc F (L- rQ) sin(2fc f (£-r . t )) 

£>9 7T ; r i>9 7T ; 

' ' - (L - r;)'" 

S 3 



3 (L~r t )V2 y ' 



(L-n, 




cos(4:k F (L 


-n)) 


(L-n 


)V2 


, sm(4:k F (L 


-n)) 



1 D. J. Derro, E. W. Hudson, K. M. Lang, S. H. Pan, J. C. 
Davis, J. T. Markert, and A. L. de Lozanne, Phys. Rev. 
Lett. 88, 097002 (2002). 

2 M. Maki, T. Nishizaki, K. Shibata, and N. Kobayashi, 
Phys. Rev. B 65, 140511 (2002). 

3 H. A. Mook, P. Dai, K. Salama, D. Lee, F. Dogan, G. 
Aeppli, A. T. Boothroyd, and M. E. Mostoller, Phys. Rev. 
Lett. 77, 370 (1996). 

4 W. E. Pickett, R. E. Cohen, and H. Krakauer, Phys. Rev. 
B 42, 8764 (1990). 

5 D. H. Lu, D. L. Feng, N. P. Armitage, K. M. Shen, A. 
Damascelli, C. Kim, F. Ronning, Z.-X. Shen, D. A. Bonn, 
R. Liang, W. N. Hardy, A. I. Rykov, and S. Tajima, Phys. 
Rev. Lett. 86, 4370 (2001). 

6 H. L. Edwards, A. L. Barr, J. T. Markert, and A. L. de 
Lozanne, Phys. Rev. Lett. 73, 1154 (1994). 

7 Y. Ando, K. Segawa, S. Komiya, and A. N. Lavrov, Phys. 
Rev. Lett. 88, 137005 (2002). 

8 D. N. Basov, R. Liang, D. A. Bonn, W. N. Hardy, B. 
Dabrowski, M. Quijada, D. B. Tanner, J. P. Rice, D. M. 
Ginsberg, and T. Timusk, Phys. Rev. Lett. 74, 598 (1995). 

9 R. Gagnon, S. Pu, B. Ellman, and L. Taillefer, Phys. Rev. 
Lett. 78, 1976 (1997). 

B. Grevin, Y. Berthier, and G. Collin, Phys. Rev. Lett. 85, 
1310 (2000). 

1 Dirk K. Morr and Alexander V. Balatsky, Phys. Rev. Lett. 
90, 067005 (2003). 

2 S. R. White, I. Affleck, and D. J. Scalapino, Phys. Rev. B 
65, 165122 (2002). 

' ! A study of Hubbard chains with alternating on-site energy 



has been done by K. Penc, and F. Mila, Phys. Rev. B 50, 
11429 (1994). 

14 S. Maekawa and T. Tohyama, Rep. Prog. Phys. 64, 383 
(2001). 

15 S. R. White, Phys. Rev. Lett. 69, 2863 (1992); Phys. Rev. 
B 48, 10345 (1993). 

16 Y. ftoh, T. Machi, C. Kasai, S. Adachi, N. Watanabe, N. 
Koshizuka, and M. Murakami, Phys. Rev. B 67, 064516 
(2003), and references therein. 

17 H. J. Schulz, Physica 235 C, 217 (1994), and references 
therein. 

18 See for example, M. Fabrizio and A. O. Gogolin, Phys. 
Rev. B 51, 17827 (1995). 

19 N. Shibata, K. Ueda, T. Nishino, and C. Ishii, Phys. Rev. 
B 54, 13495 (1996). 

20 W. Brcnig, Phys. Rep. 251, 253 (1995) and references 
therein. 

21 R. J. Cava, A. W. Hewat, E. A. Hewat, B. Batlogg, M. 
Marezio, K. M. Rabe, J. J. Krajewski, W. F. Peck, Jr. and 
L. W. Rupp, Jr., Physica C 165, 419 (1990). 

22 M. Ogata, M. U. Luchini, S. Sorella, and F. F. Assaad, 
Phys. Rev. Lett. 66, 2388 (1991). 

23 A. E. Sikkema, I. Affleck, and S. R. White, Phys. Rev. 
Lett. 79, 929 (1997); O. Zachar, Phys. Rev. B 63, 205104 
(2001). 

24 T. Mizokawa, C. Kim, Z.-X. Shen, A. Ino, A. Fujimori, M. 
Goto, H. Eisaki, S. Uchida, M. Tagami, K. Yoshida, A. I. 
Rykov, Y. Siohara, K. Tomimoto, and S. Tajima, Phys. 
Rev. B 60, 12335 (1999). 



